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PACS 96 . 60 . - j - Solar physics 
PACS 97 . 10 . Jb - Stellar activity 

PACS 05.90.+m- Other topics in statistical physics, thermodynamics, and nonlinear dynamical 
systems 

Abstract. - In this paper we analyze the behavior of the daily Sunspot Number from the 
Sunspot Index Data Center (SIDC), the mean Magnetic Field strength from the National Solar 
Observatory/Kitt Peak (NSO/KP) and Total Solar Irradiance means from Virgo/SoHO, in the 
context of the q-Triplet which emerges within nonextensive statistical mechanics. Distributions 
for the mean solar Magnetic Field show two different behaviors, with a (/-Gaussian for scales 
of 1 to 16 days and a Gaussian for scales longer than 32 days. The latter corresponds to an 
equilibrium state. Distributions for Total Solar Irradiance also show two different behaviors (ap- 
proximately Gaussian) for scales of 128 days and longer, consistent with statistical equilibrium 
and q-Gaussian for scales < 128 days. Distributions for the Sunspot Number show a g-Gaussian 
independent of timescales, consistent with a nonequilibrium state. The values obtained ("q- 
Triplet"={q s t a t,g 3en ,g re ;}) demonstrate that the Gaussian or q-Gaussian behavior of the afore- 
mentioned data depends significantly on timescales. These results point to strong multifractal 
behavior of the dataset analyzed, with the multifractal level decreasing from Sunspot Number to 
Total Solar Irradiance. In addition, we found a numerically satisfied dual relation between q s tat 
and qse-n- 
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Introduction. — One of the most important char- 
acteristics of solar magnetic activity is its periodic and 
quasiperiodic variations at different timescales. These 
cyclic phenomena include the sunspot cycle with a mean 
period of 11 years and 152-158 day periodicity observed in 
the occurrence of high-energy solar flares. Explaining their 
root-cause has become a puzzling question in solar physics. 
However, studies on the cyclic behavior of the Sun's large 
scale magnetic field have provided important information 
about different solar physical processes. Oliver et al. [l] 
established that 158-day periodicity in sunspots and high- 
energy solar flares are related to a periodic emergence of 
magnetic flux that appears near the maxima of some solar 
cycles. More recently, Knaack et al. 12] provided evidence 
for a variety of quasi-periodic oscillations in the solar pho- 
tospheric magnetic field ranging from 24 days to 22 years. 

The statistical properties of solar magnetic activity have 
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been investigated by different authors in the context of 
chaos theory [3j, multifractal analysis [2], and nonexten- 
sive statistical mechanics [5j|6]. The Tsallis' nonexten- 
sive theory is a generalization of Boltzmann-Gibbs statis- 
tical mechanics (B-G statistics) for out of thermal equi- 
librium systems, described by the entropic index q (q = 1 
in the B-G statistics; for further details: [?}{9])- It has 
been successfully applied to many complex physical sys- 
tems as solar magnetic activity [5][6], including different 
astrophysical data presenting nonlinear and nonstationary 
behavior and nonequilibrium states. Tsallis 10 recently 
proposed the existence of a three-index set (q s tat,Qsen,<lrei) , 
also known as q- Triplet, characterized by metastable states 
in nonequilibrium, where q s tat >1, Qsen <1 an d q re i >1. 
When (q sta t, Qsen, 1rei) = {^^^), the set denotes the B-G 
thermal equilibrium state. Burlaga and Vinas [5] used a 
generalized Tsallis probability distribution function to de- 
scribe the multiscale structure of fluctuations in solar wind 
speed and the magnetic field during the declining phase 
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Fig. 1: Time series of daily values for (a) Sunspot Number computed from the SIDC database, (b) Magnetic Field computed 
from NSO/KP and for (c) Total Solar Irradiance as observed by Virgo/SoHO. 



of solar cycle 23. More recently, these authors presented 
the first physical corroboration of the q-Triplet, from the 
analyses of the behavior of two sets of daily magnetic field 
strength observed by Voyager 1 in the solar wind in 1989 
and 2002 
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The magnetic activity on the solar atmosphere is by 
far more complex than such a simplified description, and 
a variety of dynamical phenomena are expected, such as 
fluctuations of the magnetic field and variability of the so- 
lar indices |12| . This solar magnetic activity is a driven 
nonlinear non-equilibrium system. These properties can 
represent a hierarchical structure in phase space, in con- 
trast to the expected behavior for B-G statistical mechan- 
ics (i.e., uniformly occupied phase space). 

The main aim of this study is to perform an analy- 
sis of the behavior of physical parameters directly reflect- 
ing the solar magnetic activity in the context Tsallis q— 
Triplet's formalism, and to compare the properties of this 
q-Triplet with those expected for a metastable or quasi- 
stationary dynamical system described by nonextensive 
statistics. The physical parameters analyzed are the solar 
indices: (i) Sunspot Number, (ii) Magnetic Field strength 
and (iii) Total Solar Irradiance, measured from observa- 
tions carried out between 1996 and 2001. These three 
datasets offer a unique possibility to analyze the relation- 
ship between the number of solar indices (e.g., sunspots 
and magnetic flux) and the noncquilibrium level, the de- 
gree of complexity, as well as the relaxation process, in 
solar activity diagnostics. The present analysis offers also 
a possibility for a check of the validity of g-Central Limit 
Theorem, the so-called g-CLT, recently conjectured by 
Umarov, Tsallis and Gell-Mann 13 , in the context of 
physical sub-systems as those associated to the solar mag- 
netic activity. 

Data and Methods. — Three diagnostics of solar 
magnetic activity were selected for this study: (i) daily 
international Sunspot Number (SN) from the Sunspot 
Index Data Center (SIDC); (ii) daily full-disk measure- 
ments of mean line-of-sight component of unsigned Mag- 
netic Field strength (MF) from the National Solar Ob- 
servatory /Kitt Peak (NSO/KP). These data do not cover 
the overall period (^20% are missing). These temporal 



gaps are reconstructed by linear interpolation of the daily 
values; and, (iii) daily means of Total Solar Irradiance 
(TSI) from Virgo/SoHO data. Observations cover the pe- 
riod January 1996 (cycle minimum) to September 2001 
(cycle maximum). This time interval corresponds to the 
increasing phase of cycle 23, which reached its maximum 
around 2001, thereafter declining until December 2007. At 
present, the Sun is known to be in solar activity Cycle 24. 
According to De Toma et al. 14 cycle 23 is magnetically 
weaker and simpler than its immediate predecessors, 21 
and 22. Figure [T] shows the variability in three solar activ- 
ity indices in this phase of cycle 23, SN and MF (without 
interpolation) and TSI. 

Virgo/SoHO data reveal the intrinsic variability of the 
Sun: from a few seconds up to the 11-year cycle. Lanza 
et al. [15] pointed out that different contributions affect 
variations of TSI: (i) shorter timescales (<2 days), as- 
sociated to active regions (solar micro- variability) ; (ii) 
longer timescales associated to active region evolution up 
to about 30 days (rotational modulation period), up to 
periodicities of ^200 days; (iii) timescales longer than 200 
days; and, (iv) very long timescales associated to an 11- 
year activity cycle. We are interested in scales between 1 
and 200 days. NSO/KP data are average values over the 
full solar disk, where several solar indices can contribute to 
its variation, such as, sunspots, faculae or plages, whereas 
SIDC data are a measure of the general state of solar mag- 
netic activity. According to De Toma et al. 14 , TSI vari- 



ability during the increasing phase of cycle 23 may be due 
to magnetic flux emergence in these indices. 

Determination of the q-Triplet. For time series S(t), 
fluctuations of increments due to its variability over 
timescale r is given as dSn(t) = S(t + r) — S(t), where r 
is the scale parameter defined as 2™ that determines the 
scale of fluctuations represented by dSn(t) and n —0, 1, 
2, ... and so on. This method has been used to describe 
magnetic field fluctuations due to solar wind on a large 
scale range [6 11 



Burlaga and Vihas [5] pointed out two factors to con- 
sider when analyzing a complex multiscale system such as 
the solar magnetic activity on from 1 hour to 1 year: (i) 
mechanical properties that differ at different scales and (ii) 
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Fig. 2: Distributions of increments for (a) SIDC, (b) NSO/KP and (c) Virgo/SoHO data. The PDFs are shifted up by a factor 
of 10 each, for the sake clarify. The inset shows the timescale r—2 n (days), from 1 to 128 days for three datasets 

, where n=0,l,2,...,7. 



the statistical probability structures of fluctuation compo- 
nents in the series at the several scales. A natural method 
for analysis of these structures is statistical mechanics. 
Two main procedures are mentioned in the literature. 
One of these is the well-known B-G statistical mechanics 
and the other Tsallis statistics. For systems studied us- 
ing nonextensive statistical mechanics according to Tsallis 
statistics, energy probability density function, sensitivity 
to the initial conditions and relaxation are described by 
three entropic indices q s tat, qsen, and q re \. These are re- 
ferred to as the q- Triplet 
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The present study ap- 
plies the same procedure used by Burlaga and Vinas [IT] 
to compute the q-Triplet from Voyager 1 data. 

The values of q s tat are derived from Probability Distri- 
bution Functions (PDFs), defined as 



PDF[dSn(t)] = A g [l + (q - l)/3 q (dSn) 2 



(1) 



where the coefficients A q , j3 q denote the normalization con- 
stants and q the functions of scale r. This entropic index 
q is related to the size of the tail in the distributions Jfjj. 
Our study uses the Levenberg-Marquardt method 16 17 



to compute the PDFs with symmetric Tsallis distribution 
from Equation (JlJ. 

The values of q sen can be obtained from the multifractal 
(or singularity) spectrum /(a), where a is the singularity 
strength or Holder exponent. Spectrum f(a) is derived 
via a modified Legendre transform, through the applica- 
tion of the MF-DFA5 method [18] . This method consists 
of a multifractal characterization of nonstationary time 
series, based on a generalization of the detrended fluctua- 
tion analysis (DFA) . The <7 sen -index denotes sensitivity at 
initial conditions. For the present purposes, we used the 



expression defined by Lyra and Tsallis 19 for the relation 
between q sen and multifractality in dissipative systems: 

1/(1 - qsen) = 1/oLmin - l/a max , (2) 

where a m i n and a max denotes the roots of the best-fit. 




t (days) 



Fig. 3: Evolution of the entropic index q sta t as a function of 
scale t, derived from fits of the PDFs in Figures [2^a),[2jb) and 
|2|c) to the g-Gaussian distribution. Where q~l, consistent 
with the Gaussian distribution. 



The value of q re i, which describes a relaxation process, 
can be computed from a scale-dependent correlation coef- 
ficient defined by 



C(r) 



{{S{t i +r)-{S(t i ))}[S{ti)-{S(t i ))}) 



([S{U)-(S{U))] 2 ) 
and for Tsallis statistics 

log C(t) = a + slogr, 

where the slope s is given by 

5 = 1/(1 - <2Ve/)- 



(3) 



(4) 



(5) 



Results. — Our results are presented in three subsec- 
tions, each one associated to the properties of one of the 
g's. 
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Fig. 4: The symbols are based on measurements of the multifractal spectrum /(a) versus a: (a) SIDC, (b) NSO/KP and (c) 
Virgo/SoHO data. All data with 1 day increment. 



On the behavior of Probability Distribution Functions 
as a function of scale: the q s t a t~index. PDFs of the 
signal S(t) were calculated to quantitatively characterize 
the stochastic process observed. The signal S(t) represents 
the SN, MF and TSI data for a timescale r ranging from 
1 to 128 days. 

PDFs for the three datasets are shown in Figure [5J The 
behavior of the PDFs entropic index changes as a func- 
tion of scale is clear. Distributions for the SN shown 
in Figure [2ja) demonstrate q-Gaussian behavior indepen- 
dent of timescales, which corresponds to a nonequlibrium 
state. However, for the MF the distributions in Figure 
[2fb) show two different behaviors. A g-Gaussian is ob- 
served for scales of 1 to 16 days and a Gaussian for scales 
larger than 32 days. The latter corresponding to an equi- 
librium state. In relation to the TSI shown in Figure |2jc) , 
distributions are approximately Gaussian at scales of 128 
days and higher, consistent with statistical equilibrium, 
but are non-Gaussian with long tails and inflection points 
on scales < 128 days. 

Differences between the PDFs of the dataset analyzed 
can also be quantitatively expressed from the nonexten- 
sivity <7 atat -index values. Figure [3] illustrates the behavior 
of the q s ta(-index distribution as a timescale function, de- 
rived from fits of the PDFs in Figure |] In the NSO/KP 
data, for timescales less than approximately 16 days the 
value for g stat ~1.3 is relatively constant. However, in 
scales larger than around 32 days, q s tat^- This is con- 
sistent with Gaussian behavior and a lack of correlation. 
For the Virgo/SoHO data, g stat ~1.65, at the average, for 
scales shorter than about 32 days. This is consistent with 
the value of g~1.75 expected for nonlinear systems, where 
the random variable is a sum of strongly correlated con- 
tributions [6lfTTl[20 



In this context, nonextensivity is 
consistent in the following aspect: Virgo/SoHO data is an 
overlap of several solar indices, such as magnetic activity 
(flares, sunspots, faculae and plages) and rotational modu- 
lation 15 21-23 . At the largest scales q s t a t decreases and 
reaches a minimum value of ~1, which is also consistent 
with Gaussian behavior and lack of correlation. 

In contrast to the data from NSO/KP and Virgo/SoHO, 



the ^tat-index from SIDC data shows a relatively anoma- 
lous behavior, with at scales shorter than about 
32 days and increasing to maximum value of ~1.6 at 128 
days. In this case, we found that SN shows an out of 
equilibrium behavior at timescales larger than about 128 
days. 



24 , due to the 



As mentioned by Kiyono et al. 
central limit theorem, Gaussian PDFs at sufficiently long 
scales are generally expected, owing to the accumulation of 
uncorrelated variations. Nonetheless, Tsallis 10 pointed 
out that during the metastable or quasi-stationary state, 
the timescales of an equilibrium stationary state can be 
reached, in many cases, infinitely later. 

On the behavior of multifractal spectra f{a): the q sen - 
index. The multifractal spectra f(a) were obtained fol- 
lowing a procedure based on the MF-DFA5 method, as 
mentioned in the Section Data and Methods. The mul- 
tifractal characterization of these data is shown in Fig- 
ure [4j These spectra f{a) calculated for datasets shows a 
wide Holder exponent interval, with a m i m =0.81±0.08 and 
a moa! =1.54±0.04 for the SIDC data, a mim = 0.71±0.04 
and a max = 1.40±0.04 for the NSO/KP data and a mim 
= 0.71±0.03 and a max = 1.33±0.07 for Virgo/SoHO 
data. For the broadness of the multifractal spectrum 
Aa=a max — a m i n , we obtained 0.73, 0.69 and 0.62 for 
the daily SN, MF and TSI, respectively. This behavior 
denotes the range of fractal exponent present in the time 
series and thus gives a measure of the degree of multifrac- 
tality of complexity [25]. These results point to a strong 
multifractal behavior with the multifractal level decreas- 
ing from SIDC to Virgo/SoHO data. 

Additional parameters can be derived from the multi- 
fractal spectrum, providing information about the fractal 
properties as, for exemple, an asymmetry in the shape of 
the f(a). Different authors (e.g.: 26 ) pointed out that a 
left or right-skewed spectrum implies high or low fractal 
exponents, respectively. In this context, the asymmetry 
reveals a relative abundance of large or small fluctuations 
in the data. Figures Qa) and (b) show that the multi- 
fractal spectra are right-skewed and slightly right-skewed, 
respectively, indicating a dominance of low fractal expo- 
nents and a relative abundance of small fluctuations in 
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SIDC and NSO/KP data. On the other hand, in Fig- 
ure |4^c), the multifractal spectrum is slightly left-skewed, 
indicating a dominance of high fractal exponents and a 
relative abundance of large fluctuations in Virgo/SoHO 
data. Therefore, in this phase of solar cycle 23, the rela- 
tive abundance of small fluctuations is more predominant 
in SIDC data than in NSO/KP and Virgo/SoHO data. 

Values of g sen -index are directly related to the instabil- 
ity of the system and to the growth of the entropy. This in- 
dex was calculated using Equation (2) based on the values 
of a in the first paragraph of the present subsection. The 
results are -0.71±0.10, -0.44±0.07 and -0.52±0.10 for 
the SN, MF and TSI, respectively. The values of q sen < 1 
denote that its distribution exhibits weak chaos (i.e., zero 
Lyapunov exponents, but not integrable |10| ) in the full 
dynamical space of the system 
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On the behavior of scale-dependent correlation: the q re i~ 
index. The value of q re ;-index can be determined from 
the scale-dependent correlation coefficient C(r) according 
to the expression given by Equation |3| . In the nonexten- 
sive theory this coefficient should decay following a power 
law given by Equation Q, with increasing r, where the 
slope s is given by Equation In general, one expects 
C(t) to evolve with the number of solar indices on the 
data and to be a function of the solar cycle. The val- 
ues of q rel obtained for SIDC, NSO/KP and Virgo/SoHO 
data are 6.9±4.81 {R 2 = 0.96), 12.11±3.98 {R 2 = 0.97) 
and 7.7±2.63 (R 2 = 0.85), respectively. The q re i values 
obtained show that the present set of data is not consis- 
tent with a power law, owing to high q re i values, particu- 
larly in the NSO /KP data. Moyano [27] suggests that the 
above procedure to calculate q re i should only be used to 
describe stochastic processes with linear correlations. In 
other words, the correlation coefficient C(r) is not a good 
alternative to conveniently describe the non-linearity of a 
sample. 

In B-G statistics, in contrast to the nonextensive theory, 
the coefficient C(r) should decrease exponentially with 
an increasing r, following a C(r) cx exp(— r/r c ) relation, 
with t c corresponding to the correlation or relaxation time. 
Considering our data as consistent with this exponential 
behavior, we found, for the parameter r c , 3.39±0.05 days 
(i? 2 = 0.90), 5.05±0.38 days (R 2 = 0.98) and 3.18±0.37 
days (R 2 — 0.97), respectively for the three datasets. A 
correlation time r c between 3 and 5 days for these datasets 
corresponds approximately to the half timescale of the re- 



laxation process obtained by Nandy et al. 28 for the Sun, 



in the magnetic fields of flare-productive solar active re- 
gions. So, if the equilibrium hypothesis of the Sun is con- 
sidered valid, then it makes sense to assume q re i — 1 for 
the three sets of data. But, is the Sun in thermal equilib- 
rium? As mentioned by Shibahashi et al. 29 , if the Sun is 



not in thermal equilibrium it takes about 10 7 years for it 
to recover its equilibrium state. Therefore, the condition 
of thermal equilibrium is a good approximation for q re i. 



Discussion and Conclusions. — We have carried 
out a new approach on the nonextensivity properties of 
solar magnetic activity from 1996-2001. The study is 
based on daily measurements of Sunspot Numbers from 
the the Sunspot Index Data Center, daily full-disk mea- 
surements of mean line-of-sight component of unsigned 
Magnetic Field strength, from the National Solar Observa- 
tory/Kitt Peak, and daily means of Total Solar Irradiance 
from Virgo/SoHO. 

The PDFs were calculated for the three datasets and 
the obtained results show that the entropic indices change 
as a function of scale. The distributions for the Mag- 
netic Field strength show two different behaviors, with a 
g-Gaussian for scales of 1 to 16 days and a Gaussian for 
scales larger than 32 days. The latter is in equilibrium 
state. The Total Solar Irradiance also show two differ- 
ent behaviors, namely, approximately Gaussian at scales 
of 128 days and larger, consistent with statistical equi- 
librium and g-Gaussian for scales < 128 days. In con- 
trast, the distributions for the Sunspot Number show a 
g-Gaussian independent of timescales and consistent with 
a nonequilibrium state for all the timescales. These results 
also confirm that the Gaussian or g-Gaussian behavior of 
the data depends significantly on the solar activity indices. 

Our results reveal that during the increasing phase 
of solar cycle 23, the multifractal spectra of the ana- 
lyzed datasets differ significantly. In particular, we found 
that SIDC and NSO/KP data have higher complexity 
than Virgo/SoHO data. In contrast, Virgo/SoHO data 
present a dominance of higher fractal exponents and a rel- 
ative abundance of larger fluctuations than for SIDC and 
NSO/KP data. 

The results above suggest that SIDC data are out 
of equilibrium more than those from NSO/KP and 
Virgo/SoHO. Again, the multifractal spectra of SIDC are 
more asymmetric, more weakly chaotic (sensitivity to the 
initial conditions) and also of higher complexity than the 
NSO/KP and Virgo/SoHO data. In fact, according to De 
Toma et al. [14] , the fluctuations in SIDC data strongly 
influence solar magnetism and TSI variability. In this con- 
text, we suggest that these properties of SN could be re- 
sponsible for the dynamics that drive the properties of the 
g-Triplet. 

Since the deviation of the three-index set 
{qstat,qsen-><lrei) from unity is a measure of the departure 
from thermodynamic equilibrium, we also calculated the 
g's indices for the present dataset, in different timescales. 
For example, for r=l day we found that the three 
solar activity indices, SN, MF and TSI, may be essen- 
tially described by the g-Triplet set (q s tat,q S en,qrei) = 
(1.31±0.07,-0.71±0.10,1), (1.21±0.06 ,-0.44±0.07,l) and 
(1.54±0.03,-0.52±0.10,1) respectively. Identical values 
of the g-Triplet (within the uncertainties) are obtained 
when we take into consideration the same dataset across 
the whole solar cycle 23. In addition, the multifractal 
spectra /(a) of the three datasets observed across the 
whole solar cycle 23 are very similar to those observed 
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during the increasing phase of solar cycle 23. This result 
suggests that solar magnetic activity tends to be in a 
quasi-stationary and metastable state for a same solar 
index. Nevertheless, when we take into account the three 
solar indices, measured simultaneously in the same period 
of the solar cycle, we find that solar magnetic activity is 
non-stationary. Such a comparison, between the whole 
cycle and its increasing phase, points in fact for a relation 
between the timescale t and the multifractal spectrum 
asymmetry, in the sense that the larger the timescale 
required to reach the equilibrium state is, the more the 
multifractal spectrum is asymmetric (see Figures [2] and 
[4j. In other words, there is a clear relation between 
indices q s tat and q sen . In short, one sees that the dual 



relation 10 30 



Qstat 



(0) 



is approximately satisfied (within the uncertainties) in the 
scenario of solar magnetic activity. This is the first time, 
to our knowledge, that a possible relationship between 
the q s tat and q sen is obtained in the study of solar phe- 
nomena. This duality reveals that the subsystems here 
analyzed (represented by the three solar indices Sunspot 
Number, Magnetic Field and Solar Total Irradiance) are 
not asymptotically scale-invariant, and that, consequently, 
they are strongly (or globally) correlated (31 . These facts 
strongly suggest that important correlations exist between 
the random variables involved in the relevant physical pro- 
cess controlling the solar activity, which still eludes us. Fi- 
nally, it is important to underline that these results are in 
clear agreement with the g-CLT 13 32 . Indeed, accord- 



ing this theorem, the q-Triplet seems to be related with 
duality relations. 
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